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Abstract 
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. The intensity noise spectrum of the hght passed through an optical microcavity is 

O calculated with allowance for thermal fluctuations of its thickness. The spectrum thus 



obtained reveals a peak at the frequency of acoustic mode localized inside the microcavity 



1^ and depends on the size of the illuminated area. The estimates of the noise magnitude 

'~~' show that it can be detected using the up-to-date noise spectroscopy technique. 

> 

On 

^ 1 Introduction 

^, 

(^ The recent progress in the area of digital high-speed spectrum analyzers of electric signals gave 

,__! rise to a strong increase of the interest to the optical noise spectroscopy [H El El HI El [6] . In the 

>■ first experimental observation of EPR of sodium atoms in the Faraday rotation noise spectrum 

/\ [Zjj a high polarimetric sensitivity [8] made it possible to detect the above spectrum using 

Cu a conventional lock-in amplifier. Application of the up-to-date digital Fourier analyzers with 

accumulating systems of data acquisition allowed one to detect magnetic resonance spectra in 

a number of solid-state objects |3j and to widen the spectral range of the noise signal detection 

up to several GHz. In this note, we describe the effect of the noise modulation of intensity of 

quasi-monochromatic light passed through an optical microcavity. The modulation arises due 

to thermal fluctuations of thickness (and, hence, of the resonance frequency) of the microcavity. 



In accordance with the calculations presented in this work, experimental observation of this 
effect proves to be possible with the aid of the digital spectrum analyzer mentioned above. In 
this paper, we primarily consider possibility of detection of this effect and only briefly discuss 
its informative potentialities. 

Let us clarify in more detail the idea of the proposed effect. Consider an optical cavity 
(Fabry-Perot interferometer) - two mirrors separated by a gap L. When the reflectivity of the 
mirrors is close to unity, the frequency dependence of the transmission coefficient ///q of such 
a cavity in the vicinity of the principal mode can be written in the form 

1 = ^ (1) 

Here, /o(-^) is the intensity of the monochromatic plane wave with the frequency u at the 
entrance (at the exit) of the interferometer and ujq = nc/L is the resonance frequency of the 
interferometer (c is the speed of light in the medium between the mirrors). The width of the 
transmission spectrum of the interferometer A is determined by the reflectivities of the mirrors 
and, for real microcavities, the finesse Q = ujq/A can reach 1000. A change in the cavity 
thickness L ^ L + ^ gives rise to a change in its resonance frequency wq — )■ cjq + Suq, which, 
at C,/L « 1, is described by the relationship \6uo/uo\ = \C,/L\. Bearing this in mind, one 
can easily make sure, using ([I]), that the relative change of the cavity thickness \^/L\ ~ l/Q, 
results in changes of its transmission coefficient of about unity and can be easily detected. In 
this case, the absolute changes of the cavity thickness ^, for typical values of parameters of the 
optical microcavities, L = 0.25/im, Q = 1000, are of the order of atomic dimensions ^ ~ 0.25 
nm. So high sensitivity of transmission of the optical cavity to variations of its parameters is 
used in physical experiments p] and allows one to set a question about possibility to detect 
the cavity transmission noise related to the cavity thickness thermal fluctuations. To evaluate 
possibility of observation of this noise, its magnitude should be compared with that of the shot 
noise of the used light source. Corresponding calculations are presented in the next section. 



2 Model calculations 

Consider an optical cavity formed by a thin layer of a medium, with the thickness L, sandwiched 
between two thin reflecting fllms. Let a monochromatic light beam with the intensity Iq and 
frequency u be incident on this cavity. Denote the area of the light spot on the cavity as D^ 



and the cavity finesse as Q- [j The coordinate system is chosen so that the plane of the cavity 
coincides with the xy plane. Then, the cavity thickness will be a function of x and y, which 
can be presented as a sum of the constant mean thickness L and small thermal fluctuation 
(^(t, X, y). We assume that the resonance frequency of the cavity is determined by its thickness 
averaged over the area D^ of the light spot, [j Now, the fluctuation 51 of the transmitted light 
intensity can be written as 



/J 
dxdy ^{t,x,y), where G{u) = — 



d f A^ \ dujQ 



TTC 

' ^0 = 7—7- 



D^duJo\A^ + {uj-iOo)y d^ 

(2) 

The greatest value of the factor G is attained at w = wq — A: G{ujq + A) = IoQ/2LD^. We 
will be interested in the intensity noise spectrum ^(z/) of the light transmitted by the cavity. 
The function ^(z/) is connected with the correlation function {61{0)61{t)) by the relationship: 

S{u)= fdte"''{SI{0)SI{t)) (3) 

Using ([2]), we obtain for {5I{0)SI{t)) the following expression: 



{6IiO)5Iit))=G\uj) / dxdydx'dy'{aO,x,y)at,x\y')) (4) 

To calculate the correlation function {^(0,x,y)^{t,x',y')) entering Eq.Q we have to: 
(i) specify a model of motion of the material of the cavity, 
(ii) obtain the appropriate Hamiltonian, and 

(iii) fulfill averaging in Eq.dll) with the thermodynamically equilibrium distribution function 

- exp[-H/kT]. 

We will describe dynamics of the cavity material in terms of small acoustic waves. We 
assume that the cavity occupies the region z E [0, L], x G [0, a], and y G [0, b]. Since the optical 
transmission spectrum depends on the cavity thickness along the z-axis, we will be interested 
only in ^-projection of displacement of the cavity material. According to the adopted model, 
this displacement is described by the acoustic field u{x,y, z) satisfying the wave equation 



dt^ 



v'Au, (5) 



^The finesse measured by the half-width of the transmission spectrum may depend on tlie transverse beam 
size D, but, at normal incidence, this dependence is weak 

^In other words, this is the frequency at which the transmission coefficient of the cavity is the greatest. 



Where f - is the velocity of longitudinal sound in the cavity material. The energy of the 
acoustic field u{x, y, z) is given by the formula 



E = dxdydzl — — h 



[P^ 7 



duX"^ /du\^ /9ti\2" 
dx) \dy J \dzj 



(6) 



Here, p - is the density of the cavity material and 7 - is the constant that described the elastic 
strain energy density. Connection of this constant with the velocity of sound v will be given 
below. The direct substitution shows that expansion of solution of Eq. ([s]) in terms of normal 
modes, satisfying the boundary conditions 



du 
dx 



0, 



x=0,a 



du 
dy 



du 

y=Ofi ClZ 



0, 



z=0,L 



(7) 



which correspond to zero strains at the bounds of the cavity, has the form: 

v^ /7rpx\ /Timy\ fiinzX 

u[t, X, y,z) = ^ Upmn COS ( I COS ( — r — I COS ( — — I , p, 11,171 mtcgcrs p,in,n > (8) 



pmn 



where the degrees of freedom of the acoustic field Upmn meet the following equations of motion 



Upmn ^pmn^pn^n^ WUCrC ^pmn ^ 



t) + (x) + (t) 



Using Eq. ([s]) , we can express the energy E ^ through the degrees of freedom 



u 



pmn- 



16 "^ 

^^ pmn 



Pulmn + 7^^ 



pmn 



TTp\^ /7[in\^ /TTn\ 



2n 



(9) 



(10) 



V 
16 



E 

pmn 



■ 2 1^ l^pmn 2 

P'^nmn ~l "75 ^ 



pmn 



To obtain the Hamiltonian corresponding to energy (10) , we have to introduce the generalized 



momenta Ppmn, conjugated with the generalized coordinates Upmn so that the equations of 
motion (|9|) acquire the form of the Hamilton equations: 



dH 



dH 



du. 



pmn 



If we set 



Ppmn 



'Ppmn 



"'pmnP '' 



dp. 



u. 



pmn 



(11) 



pmn 



and 7 = pv^ 



and express energy (10) through Ppmn and u. 



^pmn- 



H=Y. 



pmn 



4p 



pmn 



M 



+ 



Mu'L.n 



pmn '^pmn 



16 



M = pV - resonator mass, Ppmn 



M 



-u. 



pmni 



(12) 



then we can easily make sure that equations of motion ^ are equivalent to Eq. (11 ). Thus, 



Eq.(12) is the sought Hamiltonian. 



Using Hamiltonian (12), we can write the following expression for the distribution function 
o'iPpmnyUpmn} of the generalized coordinates and momenta in the state of thermodynamic 
equilibrium with the inverse temperature (3 = 1/kT: 



^XPpmny '^pmnj ^ CXp pn \Ppmni "'pmni 



(13) 



where Z is the normalizing constant. Now, if we express, using Eq. (pj), the changes of the 
cavity thickness C,(t,x,y) through the degrees of freedom Upmn 

^{t,x,y) = u{t,x,y,0) -u{t,x,y,L) = 2 }_^ Upm,2n-i cos ( 1 cos( — — 1, p,m,n> 0, 

(14) 
then, for the correlation function entering Eq. (|4]), we obtain the expression 

{aO,x,y)at,x',y'))= (15) 

^ Y- Y- / fn\ m\ f'^PX\ (T^my\ fnp'x'\ f7rm'y'\ 

= 4 2^ 2^ {Upm,2n-l{0)Up'm',2n'-l{t))COs{ 1 COS|— — Icosl 1 COS I — 1 

Since the distribution function is factorized over the degrees of freedom Upnm (i-e., the degrees 
of freedom are independent random quantities with zero mean value), the off-diagonal mean 



values in the double sum of Eq. (15) vanish, and we have 

{aO,x,y)at,x',y'))= (16) 



4 Xl("p"^-2n-l(0)Mp^,2n-l(t)) COS f j 

pmn \ CI / 



fTTmy\ fTcpx'\ fTTmy' 

COS — : — COS COS 



\ b J \ a J \ b 



Let us calculate the correlator of the type {u{0)u{t)) (the subscripts are omitted for brevity). 
Since the degrees of freedom meet the equation of motion u = —u'^u, we can write the following 
expression for u{t): 

u{t) = m(0) cosH + — sm[ut] =^ {u{0)u{t)) = (^^(O)) cos[u;t] + (m(O)m(O))^^^^ (17) 

CO CO 



The velocity u{0) entering the last term of Eq. (17) is proportional to the corresponding 
generalized momentum, which is a random quantity independent of tt. Therefore, (m(O)'u(O)) = 
0, and we have: 

{u(0)u{t)) = (^2(0)) cosH (18) 



It follows from Eqs. (12) and (13) that, if we introduce notation a = Mu'^/lGkT, then 

/2/\\ f^ f 2 217 [a d f , 21 7 1 8^^ / \ 

{u (0)) = \ - u exp -au \du = -\ -tt- / exp[-au \du = -- = ——- (19) 

V TT J V TT oa J 2a Mw^ 



Thus, using Eq.(16), we obtain 



i^/n \^u I i\\ 32/i;r ,^ cos[u}pm,2n-it\ f7ipx\ ( T^my\ fTipx \ f Timy 

(c(0, X, y)^{t, x,y)} = ——- 2^ 2 cos cos — — cos cos — — 

M ^n ^pm,2n-l ^«^ \ h J \ a J \ b 

(20) 
Now, let us take into account that the dimensions a and b are supposed to be large, i.e., 
a,b » L. In this case, we can pass from summation over p and m to integration. By denoting 
[irp/a] = A and [vrm/fo] = 5, we obtain 



lTTi\ 



and 



(e(o,x,2/)e(t,x,i/')) 



32A;T ab 

M ^ 



POO 

J2 / dAdB 



cos[a;2n-i(^,^)^] 



cos(Ax) cos{By) cos(Ax') cos{By') 

(21) 



Since Eq. (Ill) includes the correlation function (21) averaged over the area of the light spot, it 
is convenient to introduce the function F{A, B) defined as 

F{A,B)= f dxdycos[Ax]cos[By] (22) 

Then, taking into account that 

/ cos[Qt]e"''dt = 7t[6{u - Q) + 6{u + Q)], 

and using Eqs. ([s]), (21), and (22), we obtain for the sought spectral density of the noise S{iy) 
the following equation 



S{u) = G^{lo) 



32kT ab 

M ~1T 



odd .QQ 

Y, / dAdB 



6[u-Un{A,B)] 



V 



F\A,B) 



(23) 



Here, the symbol odd shows that summation is performed over odd n. Using the presence of 
(5-function, we can perform integration over B, and eventually we have 



G\u) 32kT ^ f 



A^- 



~T) 



-1/2 



F' A 



A^- 



vmy 

T) 



(24) 



Here p = M/abL is the density of the cavity material, and the integration is performed 



over the range of values of the variable A, where the radical entering Eq. (24) is real. The 



contribution of the n-th mode is evidently nonzero only at z/ > [iinv/ L\. As was already 
mentioned, the greatest value of G{u) is achieved when the interferometer is tuned to the slope 
of the resonance: 

hQ 



Gr, 



Giuo + A) 



In this case, the noise intensity is the greatest and equals 



_ jIoQf 8kT ^ 



dA 



") -A-- 



-nn\- 

t) 



-1/2 



F^ A 



"1 -A^ 



7in\- 



(25) 



To make quantitative estimates, we assume, for simplicity, that the light beam has a square 
cross section, i.e., x E [a/2 - D/2, a/2 + D/2] and y E [b/2 - D/2, 6/2 + D/2]. In this case, the 



function F(A,B) (22) can be obtained in the explicit form: 

4 



F{A,B) 



/Aa\ /AD 
— cos (^-j sin (^—, cos 



Bb\ fBD 

— sin 

2 ) \ 2 



Note that the rapidly oscillating terms (of the type cos^[i?6/2]) arising in the formula for 
F^(A, B) can be replaced by their mean values (i.e., by 1/2). With allowance for these remarks, 
we may accept that 



F2(A,S) 



sm 



AD 



sm 



BD\ 



(26) 



A2B2"" V 2 
All the above calculations refer to infinite phonon lifetime in the cavity material. This is the 



reason why the correlation function (18) does not decay. If we take into account the decay of 



the acoustic mode free vibrations and temporal evenness of the correlation function, we obtain, 
for the correlator, the following expression: 

t 



{u{0)u{t)) = (M^(0))cos[wt]exp 



T 



(27) 



where r - is the phonon lifetime. The noise spectrum S'^(z/), in this case, is obtained as a 
convolution of Eq.(25) and Lorentzian with the width R equal to the inverse phonon lifetime 



T" 



-1. 



SrM 



S{u' 



1 + [z/'r] 



du' 



(28) 



3 Possibility of observation of the microcavity noise 



One of the most popular optical microcavities is the Bragg cavity, namely, the Fabry-Perot 
interferometer comprised of two Bragg mirrors separated by a half- wave gap. In spite of the fact 



■^As we will see below, the noise spectrum appears to be localized in a rather narrow spectral region, where 
the frequency dependence of t may be neglected 



that these microcavities represent multi-layer structures, the intensity noise of the transmitted 



light can be estimated using Eqs. (25) and (28) for the following reason. It is essential for the 



single-layer model considered above that both the optical and the acoustic modes are localized 
in the same layer of the material which may serve as a resonator both for the optical and for 
the acoustic waves. A similar situation may take place in real Bragg cavities, since the optical 
Bragg mirror also has an acoustic stop-band and can efficiently reflect acoustic waves with the 
appropriate frequencies. In this case, the half-wave gap (for optical waves) between the Bragg 
mirrors can form an acoustic resonator, whose properties can be approximately described by 
the single-layer model described in the previous section. Our quantitative estimates show that, 
for a typical half-wave Bragg cavity with Aq = 2ttc/uo = 800 nm, comprised of the Ti02 and 
Si02 layers, the frequency of the lowest acoustic mode equals 10 GHz and matches the acoustic 
stop-band of the Bragg mirrors. The possibility of manufacturing of efficient acoustic Bragg 
resonators was demonstrated in [TOl [TT|, [121 IS] • 

Bearing in mind all the aforesaid, let us accept the following values of the parameters 



entering Eqs. (25) and (28): L = 0.276 /xm, Jq = 0.1 W, Q = 1000, v = 5570 m/s (SiOs), 



p = 2000 kg/m (Si02)- For the above values of the parameters, the frequency ui of the lowest 



acoustic mode {n = 1) can be estimated as vi = v/2L ^ 10 GHz. As seen from Eqs. (25) and 



(28), the noise power increases with decreasing size of the light spot D. In our calculation, we 



accept that D = 10 /im. This size of the spot does not contradict the above finesse Q = 1000. £j 



To estimate the decay time r of acoustic vibrations, entering Eq. (28), we can take into account 
that the finesse of the quartz resonators at frequencies of around 10* Hz may reach 10^ — 10^. 
For the frequencies ~ lO^'' Hz, we are interested in, the finesse of the acoustic vibrations is 
expected to be lower. Thus, for our estimates, we will accept the finesse of the acoustic mode 
to be ~2000. In this case, the quantity r is determined by the relationship uit = 2000. The 
possibility of manufacturing of acoustic microcavities with a finesse of around 1000 in the THz 
range was reported in [14] . 

To estimate possibility of observation of the intensity noise of the light transmitted by the 



microcavity, the magnitude of this noise (determined by Eqs (25) and (28)) should be compared 
with the shot noise of the light source 

Ssn = lo^OJo (29) 



^Futher decrease of D by light focusing can be accompanied by a decrease of the finesse, which occurs due 
to increasing uncertainty in the angle of incidence. However, for normal incidence of the focused beam, this 
effect is not strong. 



The results of this comparison are shown in Fig. 1. This figure shows the calculated noise 
spectra in the region of the lowest acoustic mode (n = 1) for the infinite (oscillating curve) and 
finite (smooth curve) phonon lifetime. Horizontal line shows the level of shot noise of the used 



light Eq. (29). Since observation of the light shot noise with the aid of the up-to-date spectrum 
analyzers does not encounter any problems, then, as seen from Fig. 1, we have all the grounds 
to believe that the above effect can be detected even if the results of our calculations appear 
to be overestimated by 1-2 orders of magnitude, [j 



4 Conclusions 

In this paper, we presented calculations of spectral power density of the light transmitted by 
a microcavity. It is shown that this power density reveals a peak at the frequency of acoustic 
vibrations of the cavity - the effect similar to Raman scattering. The quantitative estimates 
are made which show that the noise arising due to the mechanism considered in the paper can 
be detected using the up-to-date noise-spectroscopy technique. 

Without entering into details of possible informative capabilities of this technic, note only 
that thermal vibrations are usually considered as a spurious factor that restricts operational 
stability of devices (see, i.g.,[T5]. The above calculations show, however, that the intensity noise 
spectrum of the light transmitted by the cavity contains information about the spectrum of 
acoustic vibrations of the structure. Measuring variation of this spectrum versus the light spot 
diameter will allow one to judge about plausibility of the used simple model which neglects, 
in particular, disorder of the real structure and possible localization of the acoustic waves. 
Observation of the correlation function of the noise for two spaced light beams will possibly 
allow one to evaluate the localization radius of acoustic vibrations of the structure. 
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^Spectral region of the modern digital spectrum analyzers is restricted to the frequencies of around 1-2 GHz. 
For this reason, for detection of the noise signal with the frequencies 10 GHz, considered in this paper, one 
has to perform the appropriate transfer of the spectrum. A similar task is solved in the systems of satellite TV 
with the aid of heterodyne converters, which can also be used in this case. 
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Figure 1: The intensity noise spectrum of the hght transmitted by the cavity. The values of 
the parameters are: thickness of the microcavity L = 0.276 /xm, finesse of the transmission 
spectrum Q = 1000, and the hght beam intensity Jq = 0.1 W. The calculation is performed for 
the spectral region of the lowest acoustic mode of the cavity uq ~ 10 GHz. Oscillating curve 
5'(z/) corresponds to the noise spectrum at infinite phonon lifetime. Smooth curve Sr{iy) - the 
same at 27rz/oT = 2000. Horizontal line shows the shot noise level of the light. 
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